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Preface 


This  thesis  is  one  of  a  continuing  series  of  studies 
conducted  at  the  Air  Force  Institute  of  Technology  (AFIT) 
to  increase  the  understanding  of  frequency  response 
characteristics  of  pneumatic  transmission  lines.  Better 
models  are  sought  to  enhance  the  prediction  capabilities 
for  various  flow  conditions  and  line  configurations.  The 
purpose  of  this  study  is  to  investigate  the  effects  of 
turbulent  mean  flow  and  line  terminations  on  the  dynamic 
response  of  circular  pneumatic  lines. 

Previous  studies  performed  at  AFIT,  notably  those  of 
Moore,  Vining,  Malanowski  and  Briski,  provided  invaluable 
information  and  guidance  pertaining  to  the  experimental 
and  analytical  aspects  of  this  thesis. 

I  wish  to  thank  a  number  of  individuals  whose  support 
made  this  study  possible.  No  matter  how  busy.  Dr.  M.  E. 
Franke  was  always  available  for  consultation  and  technical 
advise.  Dr.  W.  Elrod  and  Capt.  W.  Cox  contributed 
welcomed  analytical  and  experimental  assistance.  A 
special  thanks  to  my  wife  Cindy  who  managed  the  household 
and  keep  our  family  together  through  some  difficult 
periods . 


Phillip  G.  Wilkins 
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Symbol  Description 


Units 


A 

line  cross-sectional  area 

in 

a 

adiabatic  speed  of  sound 

in/ sec 

B 

phase  angle 

radians 

C 

capacitance/unit  length 

in'*/  lbf 

Ca 

adiabatic  capacitance/unit  length 

in4/lbf 

Ci 

isothermal  capacitance/unit  length 

in4/ lbf 

CT 

turbulent  capacitance/unit  length 

in4/ lbf 

D 

inside  line  diameter 

in 

f 

friction  factor 

dimensionless 

G 

conductance/unit  length 

in4/ lbf-sec 

g 

pressure  gain 

dimensionless 

gdb 

pressure  gain  in  decibels 

decibels 

ht 

temperature  distribution  factor 

dimensionless 

hv 

velocity  distribution  factor 

dimensionless 

j 

m 

• 

o 

r- 1 

i 

dimensionless 

k 

ratio  of  specific  heats 

dimensionless 

klt 

inertance  constant 

dimensionless 

L 

inertance/unit  length 

lbf-sec^/ ink 

La 

adiabatic  inertance/unit 
length 

lbf-sec^ / in^ 

L>p 

turbulent  inertance/unit 
length 

lbf-sec^/ ink 

M 

mass  flow  rate 

lbm/ sec 

n 

polytropic  exponent 

dimensionless 

Symbol  Description 


Units 


pressure 

ambient  pressure 
inlet  pressure 
pressure  drop  across  load 
receiving  pressure 
sending  pressure 
static  line  pressure 
volumetric  flow  rate 
resistance/unit  length 
AC  resistance/unit  length 
DC  resistance/unit  length 
gas  constant  for  air 
Reynolds  number 
reflection  coefficient 
inside  line  radius 
ambient  temperature 
velocity 
volume 

distance  in  downstream  direction 

admittance/unit  length 

impedance/unit  length 

characteristic  impedance/unit 
length 

input  impedance 
load  impedance 
radiation  impedance 


psia 


psia 


ps  la 


psig 


psig 


psig 


lbf-sec/ in^ 
lbf-sec/ ink 

lbf-sec/ in6 

in2/sec2/R 


dimensionless 


dimensionless 


in/ sec 


in^ / lbf-sec 
lbf-sec/ ink 

lbf-sec/ in6 
lbf-sec/ in5 
lbf-sec/ in5 
lbf-sec/ in5 


v. •_>.  f, 


Symbol 

Description 

Units 

a 

attenuation/unit  length 

neper/ in 

/? 

phase  shift/unit  length 

rad/ in 

r 

propagation  operator/unit  length 

1/in 

p 

dynamic  viscosity 

lbf-sec/ in2 

V 

kinematic  viscosity 

in^/sec 

Q 

ncndimensional  frequency 

dimensionless 

Q  B 

nondimensiona 1  break  frequency 

dimensionless 

0) 

angular  frequency 

rad/ sec 

0) 

'"v 

characteristic  frequency 

rad/ sec 

P 

density 

lbf-sec^/in^ 

a 

(Prandtl  number)  *05 

dimensionless 

ea 

radiation  function,  real 

dimensionless 

* 

radiation  function,  imaginary 

dimensionless 
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Abstract 


Experimental  frequency  response  data  were  obtained  for 
straight  pneumatic  transmission  lines  of  circular 
cross-section.  Three  24  in.  long  lines  with  inside 
diameters  of  0.041,  0.119  and  0.195  in.  were  tested  at 
Reynolds  numbers  varying  from  0  to  15000. 

The  experimental  data  were  compared  with  theoretical 
predictions  obtained  from  the  distributed  parameter,  small 
signal  equations  for  the  blocked  and  laminar  mean  flow 
cases.  Several  models  were  used  to  predict  the  effects  of 
turbulent  mean  flow.  The  prediction  capability  and 
limitations  of  two  turbulent  flow  models,  that  include  the 
effects  of  the  Reynolds  number,  but  not  frequency,  were 
investigated  and  discussed. 

Methods  for  determining  the  end  impedance  of  open  and 
orifice  terminated  lines  with  mean  flow  were  studied.  The 
relatively  simple  approach  of  reflection  coefficients  gave 
satisfactory  results  for  both  choked  and  unchoked  flow  in 
orifices.  The  radiation  impedance  associated  with  open 
lines  was  found  to  be  negligible  for  all  test  cases. 


DYNAMIC  RESPONSE  OF  FLUID  TRANSMISSION  LINES 


I .  Introduction 


Background 

A  periodic  pressure  signal  propagating  in  a  fluid 
filled  line  exhibits  changes  both  in  magnitude,  defined  as 
gain,  and  in  phase  angle,  known  as  phase  shift,  along  the 
length  of  the  line.  The  signal  may  be  information,  as  in 
a  fluidics  circuit  or  instrumentation  system,  noise  in  an 
air  duct  or  even  an  undesired  pulsation  due  to  a  pump  or 
valve.  It  is  essential  to  many  system  designs  to 

accurately  predict  the  gain  and  phase  shift  to  assess  the 
impact  on  desired  operation. 

The  gain  and  phase  shift  are  dependent  on  such 

parameters  as  line  geometry,  signal  frequency,  fluid 


medium 

and 

flow 

conditions.  Early 

investigations  were 

1 imited 

to 

the 

field  of  acoustics. 

Instrumentation 

applications 

led 

to  later  studies,  many  of  which  involved 

liquids.  The  relatively  young  field  of  fluidics  brought 
renewed  interest  in  the  dynamic  response  of  fluid 
systems.  Many  studies  have  been  conducted  over  a  wide 
range  of  signal  frequencies  using  circular,  annular  and 
rectangular  lines  with,  primarily,  air  or  water  as  the 


fluid  medium. 


The  dynamic  response  of  blocked,  ,\o  mean  flow,  lines 
has  been  accurately  predicted  by  t,.e  equations  developed 
by  Nichols  (Ref  1)  and  later  modified  by  Krishnaiyer  and 


Lechner  (Ref  2) -  These  equations  also  apply  to  lines  in 
which  the  signals  are  imposed  on  a  laminar  mean  flow. 

Brown,  Margolis  and  Shaw  (Ref  3) ,  Moore  (Ref  4)  and 
Briski  (Ref  5),  among  others,  showed  that  application  of 
these  equations  to  cases  involving  turbulent  mean  flow 
lead  to  good  agreement  at  higher  frequencies,  but  under 
predicted  the  attention  at  lower  frequencies.  Brown,  et 
al.  (Ref  3)  developed  a  somewhat  complicated  procedure  for 
predicting  turbulent  flow  response  by  considering  two  and 
three-region  boundary  layer  models.  Moore  (Ref  4) 
described  three  frequency  regimes  based  on  boundary  layer 
thickness  and  used  a  constant  inertance,  resistance  and 
capacitance  (LRC)  model  to  predict  the  increased 
attenuation  of  the  turbulent  flow,  compared  to  the  laminar 
case,  at  lower  frequencies.  Briski  (Ref  5)  modified  the 
constant  LRC  model  to  account  for  the  effects  of  mean  flow 
on  the  resistance. 

All  fluid  lines  are  terminated  in  some  fashion. 
Vining  (Ref  6)  studied  means  to  determine  the  impedance  of 
orifices  with  mean  air  flow.  The  proper  determination  of 
the  terminal  impedance  was  shown  to  be  an  important  factor 
in  modeling  the  line. 


Objectives 

The  following  objectives  were  established  to  study  the 
dynamic  response  of  pneumatic  transmission  lines. 


1.  Investigate  means  of  determining  the  impedance  of 
various  termination  and  flow  combinations. 

2.  Modify  existing  computer  program  to  incorporate 
revised  theoretical  models  and  techniques  for 
determining  the  pressure  drop  and  Reynolds  number 
of  a  line. 

3.  Experimentally  determine  gain  and  phase  shift  for 
various  test  line  configurations  and  flow 
conditions . 

4.  Present  simple  models  to  predict  the  response  of 
lines  with  turbulent  mean  flow. 

5.  Assess  response  predictions  for  blocked,  laminar 
and  turbulent  mean  flow  cases  by  comparison  with 
experimental  data. 


II.  Theory 


Electrical  Analogy 

Electrical  transmission  theory  is  well  developed  and 
can  be  found  in  great  detail  in  King  (Ref  7)  among 
others.  Electrical  theory  can  be  applied  to  the  study  of 
fluid  transmission  lines  by  considering  pressure  to  be 
analogous  to  potential  and  volumetric  flow  rate  to 
represent  the  current.  Ohm's  Law  then  becomes 


4  P  =  2Q 


Q  =  y  4  P 


where  z  is  the  complex  series  impedance,  y  is  the  complex 
admittance,  4P  is  the  pressure  change  and  Q  is  the 
volumetric  flow  rate.  The  above  equations  expressed  as  a 
function  of  the  length,  x,  are 


dP/dx  =  ZQ  =  (R  +  j  uL)  Q 


dQ/dx  =  YP  =  (G  +  j  <JC )  P 


where  R,  L,  G  and  C  are  the  resistance,  inertance, 
conductance  and  capacitance  per  unit  length  of  line.  The 
analogous  circuit  for  a  fluid  transmission  line  is  shown 
in  Fig  1. 


Blocked  and  Laminar  Flow  Lines 

Values  for  Z  and  Y  may  be  obtained  for  circular 

blocked  lines  and  those  with  laminar  mean  flow  by  using 

the  equations  of  Nichols  (Ref  1) .  This  approach  involves 
the  use  of  Bessel  functions  and  becomes  rather 
complicated.  Krishnaiyer  and  Lechner  (Ref  2)  present 

modifications  that  have  been  shown  to  give  accurate 

predictions  for  Z  and  Y  for  the  frequency  range  u>  / 10 
where  uv  is  the  characteristic  frequency  defined  as 

=  32  V/D2  (5) 

where  D  is  the  line  diameter. 

The  modified  equations  are  given  as 

Z  =  (8  7T/A2)  (DR)  +  j  (  (U»p/A)  +  (87T  fj/A2)  (DL)  )  (6) 

and 

Y  =  -  [  cj(k-l)  kp]  (DG)  +  j  u\  kP  +  (k-1 )  (kp)  (DC)  1  (7) 

(DC)2  +  ( DG )  2  (DC)2  +  (DG)2 

where  the  quantities  DR,  DL,  DG  and  DC  are  described  in 
terms  of  velocity  and  temperature  parameters,  are 
respectively 

DR  =  3/8  +  h„  /4  +  3/ (8  h„  )  (8) 

DL  =  h„  /4  -  15/  ( 64h„  )  (9) 


DG  =  ht/2  -  1/ ( 4ht ) 


(10) 


DC  =  1/4  +  ht/2  +  1/ ( 4ht ) 


(ID 


where 


h„  =  2(ulu„  )0-5 


(12) 


and 


ht  =  <rhv 


(13) 


The  line  may  now  be  described  by  a  characteristic 
impedance  defined  as 


Zc  =  (Z/Y) *05  (14) 

and  a  propagation  operator 

r=  (ZY)  0*5  =  a  +  j  0  (15) 

where  a  and  are  the  attenuation  and  phase  shift  per 
unit  length.  These  parameters,  along  with  the  terminal 
impedance  and  line  length  will  be  used  later  to  determine 
gain  and  phase  shift. 

Kirshner  and  Katz  (Ref  8)  give  the  input  impedance  for 
the  general  case  of  a  line  with  a  terminal  impedance  ZL, 
as  shown  in  Fig  2,  as 

=  pin/Qin  =  ^c^L  cosh  r£  +  Z^sinh/-1^ )  /  (16) 

(ZLsinhr^  +Zccos  h  r£) 

The  transfer  function  relating  the  value  of  P  at  any  x 
along  the  line  to  the  inlet  pressure  becomes 


7 


P/Pin  =  [sinh/^(l-x/^  )  +  ( ZL/ZC)  cosh/"./  (1-x/ Jt  )  ]  /  (17) 

[sinhTji  +  (Zl/Zc)  cosh/-#  ] 

The  pressure  gain  is  the  magnitude  of  equation  (17) 

9  =  | P/Pin  1  (18) 

or  expressed  in  decibels 

gdb  =  20  LOG10|P/Pinj  (19) 

The  corresponding  phase  shift  is  the  angle  formed  in 
the  complex  plane  by  the  real  and  imaginary  parts  of  the 
transfer  function 

B  =  TAN-l[Im(P/Pin)/Re(P/Pin) ]  (20) 

Turbulent  Flow 

Equations  (6)  and  (7)  yield  accurate  predictions  for 
the  response  of  blocked  and  laminar  flow  cases.  Attempts 
to  extend  their  application  to  turbulent  mean  flow 
underestimate  the  experimentally  observed  attenuation, 
(Ref  3-5) . 

The  laminar  steady  state  fully  developed  flow  has  a 
parabolic  velocity  profile  with  a  boundary  layer  thickness 
equal  to  the  radius.  The  dynamic  boundary  layer,  created 
when  a  periodic  signal  is  imposed  on  the  steady  state,  or 
mean,  flow,  becomes  thinner,  that  is  the  gradients  near 
the  wall  increase,  as  the  frequency  increases.  Therefore, 
for  all  frequencies  of  practical  interest,  the  dynamic 


layer  is  much  thinner  than  the  steady  state  and  therefore 
dictates  the  response  behavior.  Equations  (6)  and  (7)  are 
seen  to  be  dependent  on  the  signal  frequency  but  not  the 
Reynolds  number. 

For  turbulent  flow,  however,  the  steady  state  velocity 
profile  becomes  dependent  on  the  Reynolds  number.  As  the 
Reynolds  number  increases,  the  boundary  layer  becomes 
thinner  until  in  the  limiting  case  of  an  infinite  Reynolds 
number  a  uniform  profile,  or  slug  flow,  is  obtained.  When 
a  periodic  signal  is  imposed  on  turbulent  flow,  the 
dynamic  layer  may  or  may  not  be  thinner  than  the  steady 
state  layer.  Since  the  response  is  dictated  by  the 
thinner  layer,  one  would  expect  to  observe  different 
behavior  due  to  frequency  and  Reynolds  number 
combinations . 

A  plot  of  nondimensional  attenuation  versus 
nondimensional  frequency,  Q  ,  presented  by  Brown,  et  al. 
(Ref  3)  is  shown  in  Fig  3.  Future  use  of  the  terms 
attenuation  and  frequency  will  refer  to  these 
nondimensional  parameters.  Moore  (Ref  4),  as  suggested  by 
Brown,  et  al.  (Ref  3),  divides  the  frequency  spectrum  into 
three  regions;  low,  mid  and  high. 

In  the  low  frequency  region,  the  dynamic  boundary 
layer  is  thicker  than  the  steady  state  layer.  The  high 
frequency  region  exhibits  the  same  behavior  as  the  laminar 
flow  due  to  the  fact  that  the  dynamic  layer  has  grown  much 
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thinner  than  the  steady  state  layer.  The  mid  frequency 
range  is  characterized  by  dynamic  and  steady  stare  layers 
of  comparable  thickness. 

Moore  (Ref  4)  applied  a  constant  LRC  model  to  both  the 
low  and  mid  frequency  regions.  This  approach  used  a 
steady  state  (DC)  resistance 

Rdc  =  f  Re  ft/ 2  A  D2  (21) 

which  by  using  the  Blasius  resistance  formula,  as 
presented  by  Schlichting  (Ref  9) ,  for  turbulent  flow  in 
smooth  lines 


f  =  0.3164  Re-0-25  (22) 

was  written  as 

RDC  =  0.3164  Re0-75  fl  /  (2  A  D2)  (23) 

The  capacitance  was  assumed  to  be 

C  =  (k/n) Ca  (24) 

where  n  is  the  polytropic  exponent  obtained  from  the 
thermodynamic  process  description 

PVn  =  constant  (25) 

and  Ca  is  the  adiabatic  capacitance  defined  as 

Ca  =  A/ (kP)  (26) 
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The  value  of  n  varies  from  k  at  high  frequencies  to 
1.0  for  the  isothermal  process  though  to  exist  at  low 
frequencies  (Ref  4) . 

The  adiabatic  inertance  used  is  defined  as 

La  =  pi  A  (27) 

As  the  model  name  implies,  the  conductance  was 
neglected.  This  choice,  as  well  as  that  of  isothermal 
capacitance  and  adiabatic  inertance,  was  based  on  analysis 
of  the  laminar  case  at  low  frequency  (Ref  4) 

Using  equations  (3)  and  (4) ,  the  turbulent  impedance 
and  admittance  can  be  expressed  as 

2  =  ^DC  +  j  OJ  La  ( 28 ) 

and 

Y  =  j  o>  C  (29) 

The  gain  and  phase  shift  can  now  be  determined  through 
the  use  of  equations  (14)-(28).  The  laminar  model,  as 
shown  experimentally  by  Brown,  et  al.  (Ref  3)  and  Funk 
and  Wood  (Ref  10),  among  others,  accurately  describes  the 
response  behavior  at  frequencies  above  the  break 
frequency,  Q b ,  which  is  defined  as  the  frequency  at 
which  the  attenuation  predictions  of  the  laminar, 
equations  (6)  and  (7) ,  and  the  turbulent,  in  this  case  the 
constant  LRC ,  models  are  equal.  At  frequencies  below  the 


break  frequency,  however,  the  laminar  model  under  predicts 
the  attenuation. 

As  shown  by  Briski  (Ref  5),  the  attenuation  predicted 
by  the  constant  LRC  model  is  also  lower  than  that  observed 
experimentally.  Modifications  to  the  constant  LRC  model 
will  be  presented  in  an  effort  to  improve  the  prediction 
capability  in  the  low  and  mid  frequency  regions.  Although 
L,  R,  C  and  G  have  been  shown  to  be  frequency  dependent 
(Ref  4),  this  study  will  employ  only  constant  values  to 
maintain  the  simplicity  of  the  constant  LRC  model. 

The  mid  frequency  region  of  Fig  3  is  characterized  by 
a  rather  flat  attenuation  curve  dependent  upon  the 
Reynolds  number  only.  Therefore,  only  the  Reynolds  number 
dependency  of  each  of  the  components  will  be  investigated. 

To  account  for  the  variation  in  resistance  due  to 
Reynolds  number,  a  linearized  or  incremental  (AC) 
resistance  was  proposed  (Ref  11) 

rAC  =  RDC  +  Re  H  9  / 9  Re ) Rq^ )  (30) 

which,  using  the  DC  resistance  of  equation  (23),  reduces 
to 

rAC  =  1*75  RDC  (31) 

Based  on  the  effect  of  the  velocity  profile  on  the 
inertance,  Moore  and  Fr.-'nke  (Ref  12)  give  the  inertance 
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lT  =  kLT  La 


(32) 


where  is  the  inertance  constant 

Klt  =  (f°*5+l) (f°-5+2) 2/ [4(2f°-5+l) ]  (33) 

Klt  is  seen  to  vary  from  1.33  for  laminar  flows  with 
parabolic  velocity  profiles  to  1.0  as  the  Reynolds  number 
goes  to  infinity  and  the  profile  approaches  that  of  slug 
flow  (Ref  12).  Since  the  friction  factor  is  Reynolds 
number  dependent,  the  inertance  is  also  shown  to  be  a 
function  of  the  Reynolds  number. 

The  conductance  accounts  for  the  effects  of  radial 
heat  transfer  of  a  gas.  The  increased  attenuation  due  to 
the  radial  heat  transfer  is  well  understood  for  laminar 
flow  and  is  incorporated  in  equations  (6)  and  (7).  Moore 
(Ref  4)  assumed  the  conductance  for  turbulent  flow  to  be 
negligible  based  on  the  small  values  observed  for  laminar 
flow  with  small  signal  frequencies.  Brown,  et  al.  (Ref  3) 
presents  the  conductance  for  turbulent  flow  as 

G  =  (  (k-l)CaLaw2)  /  (RAC  +  j  wLT)  (34) 

which  reduces  to  the  following  for  cases  where  the 
quantity  (R^q/io  )  is  much  smaller  than  LT 


G  =  (  (k-1)  RACCaLa/LT2)-j  (j  (  (k-l)CaKLTLa2/LT2) 


(35) 


The  imaginary  part  of  the  conductance  can  be  combined 
with  the  isothermal  capacitance  to  form  a  turbulent 
capacitance 


CT  =  kCa[l-(k-l) / (KLTk) ]  (36) 

which  is  seen  to  be  Reynolds  number  dependent  through  the 
inertance  constant-  The  real  part  of  the  conductance  is 
Reynolds  number  dependent  by  virtue  of  both  the  resistance 
and  inertance  constant. 

Two  models  for  turbulent  mean  flow  can  now  be 
presented.  The  modified  constant  LRC  model  consists  of 
the  turbulent  inertance  from  equation  (32),  AC  resistance 
of  equation  (31)  and  turbulent  capacitance  given  by 
equation  (36)  .  The  constant  LRCG  model  adds  the  real  part 
of  the  conductance  from  equation  (35)  to  the  components  of 
the  modified  constant  LRC  model.  All  components  in  both 
models  are  seen  to  be  functions  of  the  Reynolds  number  not 
of  frequency. 

Terminal  Impedance 

The  steady  state  DC  terminal  impedance  of  a  fluid  line 
can  be  expressed  as 


ZL  =  pl/<2l  +  3  0 


(37) 


where  Pl  is  the  pressure  drop  across  the  termination. 
Franke,  Malanowski  and  Martin  (Ref  13),  among  others,  have 
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shown  it  more  applicable  to  use  an  incremental  or 
linearized  AC  expression  when  dealing  with  dynamic  systems 

ZL  =  (  4PL)/  AQl  +  j  0  (38) 


For  blocked  lines,  Ql  =  0  and  Zl  =  ®  ;  while  for 

the  other  extreme  case  of  an  open  line,  Pl  and  Z^  are 
taken  to  be  zero. 

The  effect  of  the  terminal  impedance  on  the  pressure 
gain  can  be  seen  in  Table  I,  from  (Ref  8),  in  which  the 
transfer  functions  for  the  special  cases  of  blocked,  open 
and  matched  lines  are  given.  The  matched  line  has  a 
terminal  impedance  equal  to  its  characteristic  impedance 
and  the  line  appears  infinite  in  length  to  the  input; 
therefore,  no  waves  are  reflected,  no  standing  wave 
pattern  is  created  and  only  exponential  decay  results. 

Acoustic  theory  states  the  existence  of  a  finite 
impedance  associated  with  the  open  line.  The  air  in  the 
open  end  of  the  line  can  be  thought  of  as  a  massless 
piston  radiating  some  energy  into  free  space  while 

reflecting  some  back  into  the  line.  Morse  and  Ingard 
(Ref  14)  give  this  radiation  impedance  as 


ZR  =  (  p  a/A)  (  0o-j  XJ 


(39) 


where  Q  and 

o 

approximating  the 
of  a  flanged  exit, 


X  are  termed  radiation  functions.  By 
values  of  these  functions  for  the  case 
defined  as  a  flush  opening  in  an 
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infinite  wall,  the  real  portion  is  assumed  to  be 
negligible  while  the  imaginary  part  becomes 


X  =  0.4  (D  u/a) 


The  resultant  expression  for  the  terminal  impedance 
of  an  open  line  is  taken  to  be 


Zr  =  -  j  0 . 4  (  pD  cj/A) 


Kinsler  and  Frey  (Ref  15)  discuss  using  reflection 
coefficients  to  obtain  impedance  values.  Referring  to  the 
notation  of  Fig  4 


ZL  =  (  p  a/ A) i  [ (1+RC) / (1-RC) )  +  j  o 


where  RC  is  the  reflection  coefficient  for  the  sudden  area 
change  from  line  1  to  line  2.  For  the  case  where  line  2 
is  terminated  by  its  characteristic  impedance 


RC  =  (A]_  -  A2)/(A1  +  a2: 


The  use  of  reflection  coefficients  simplifies  the 
difficult  problem  of  determining  the  impedance  of  an 
orifice  with  compressible  flow.  For  cases  in  which  the 
flow  was  choked  by  the  orifice,  the  entire  orifice 
diameter  was  assumed  to  contain  sonic  flow.  Katz,  Hausner 
and  Eisenberg  (Ref  16)  indicate  that  reflected  waves  from 
a  boundary,  which  propagate  upstream  at  the  speed  of 
sound,  will  be  swept  downstream  by  the  sonic  flow. 


Figure  4  Sudden  Area  Change 


Therefore,  the  line  will  appear  infinite,  or  matched,  to 
the  input  and  use  of  equation  (43)  is  justified. 

Determination  of  the  impedance  for  the  case  where  the 
pressure  ratio  across  the  orifice  is  insufficient  to 
choke  the  flow  is  not  as  straightforward.  For  these 
cases,  the  orifice  was  considered  as  an  infinitesimally 
short  open  line,  with  a  corresponding  reflection 
coefficient  of  -1.  The  remainder  of  the  interface 
presents  a  blocked  line,  reflection  coefficient  of  +1,  to 
the  input.  The  two  values  are  then  averaged  by  area  to 
obtain  the  composite  reflection  coefficient 

RC  =  (A1-2A2)/(A1+A2)  (44) 

Equation  (43)  is  actually  obtained  in  the  same  fashion  by 
using  the  fact  that  the  reflection  coefficient  for  a 
matched  line  is  0. 

Appendix  B  outlines  various  orifice  impedance  models 
considered  for  both  choked  and  unchoked  cases. 
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calculated  by  the  laminar  and  selected  turbulent  model  for 
each  frequency.  The  attenuations  are  then  compared  and 

the  largest  is  used  to  find  the  gain  and  phase  shift.  As 
a  result,  the  turbulent  model  predicts  the  response  up  to 
the  break  frequency  at  which  point  the  program  shifts  to 
the  laminar  model. 

The  program  divides  the  line  into  a  number  of  short 
line  segments  over  which  it  is  reasonable  to  assume 
constant  average  values  of  certain  properties  to  simplify 
calculations.  Since  the  input  pressure  to  each  of  the 
segments  is  required  by  the  main  part  of  the  program, 
modifications  were  made  to  obtain  the  pressure  along  the 
line  given  a  measured  line  inlet  pressure  and  mass  flow 
rate,  which  yields  the  Reynolds  number. 

Schlichting  (Ref  9)  gives  the  pressure  drop  in  a  line 
as 

A  P  =  f  (A  x)  p  (u2/2D)  (45) 

For  laminar  flow  in  smooth  lines,  the  friction  factor 
is 


f  =  64/Re  (46) 

The  pressure  change  for  laminar  flow  then  becomes 

A  P  =  -3  2  (  A  x)  u2  (Re/  p  D3)  (47) 

Using  the  Blasius  resistance  formula,  the  pressure 
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change  for  turbulent  flow  in  smooth  lines  is  seen  to  be 

AP  *  -0 . 1582  (  A  x)  //2(Re1-75/p  D3)  (48) 

Assuming  constant  temperature  and  Reynolds  number  in 
the  short  line  segment  of  length  x,  equation  (47)  is 
integrated  to  give  for  laminar  flow 

AP  =  Pin-[Pin2+64(d  X)  p2ReTa (Rg/D3) ] °-5  (49) 

where  P-[n  is  the  inlet  pressure  to  the  line  segment  and 
Rg  is  the  gas  constant  for  air.  The  corresponding 
equation  for  turbulent  flow  is  found  from  equation  (48)  to 
be 

AP  =  Pin-(Pin2  +  0.3164  (Ax)  p2  Re1-75Ta  (Rg/D3)  ]0-5  (50) 

Using  the  measured  pressure  at  the  beginning  of  the 
line,  the  program  marches  along  the  line  and  calculates 
the  inlet  pressure  for  each  line  segment. 

A  complete  listing  of  the  program  is  presented  in 


Appendix  C. 


III. 


Experiment 


Apparatus 

Experimental  data  for  this  study  consisted  of  pressure 
measurements  taken  at  the  inlet  and  exit  planes  of  several 
stainless  steel  lines.  The  dynamic  pressures  at  these 
locations,  termed  the  sending  and  receiving  pressures, 
were  a  result  of  the  standing  wave  pattern  created  by  the 
periodic  signal  imposed  upstream  of  the  inlet  plane. 
Three  lines  were  used,  each  24  in.  long  and  having  inside 
diameters  of  0.041,  0.119  or  0.195  in.  Each  line  was 
supported  at  both  ends  by  plexiglass  blocks  which  were 
machined  to  provide  mounting  locations  for  pressure 
transducers,  air  supply  connections  and  interchangeable 
termination  devices.  Fig  5  depicts  the  test  lines  for 
which  the  dimensions  can  be  found  in  Table  II.  Special 
devices  were  used  in  each  of  the  three  lines  to  create 
blocked,  open  and  two  orifice  terminated  conditions. 
Details  of  the  devices  are  presented  in  Fig  6  with 
corresponding  dimensions  listed  in  Table  II. 

Referring  to  Fig  7,  the  sending  wave  analyzer  provided 
an  input  signal  of  known  frequency,  as  verified  by  a 
frequency  counter,  to  the  pneumatic  driver  assembly.  The 
driver  assembly  amplified  and  conditioned  the  input  and 
imposed  sinusoidal  pressure  signals  on  the  mean  air 
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Number 


imensions  in  inches 


Static  Pressure  DC  Power  Supply 


flow  passing  through  the  driver.  Piezoelectric 
transducers  were  used  to  measure  the  dynamic  sending  and 
receiving  pressures,  Ps  and  Pp.  The  transducer 
outputs  were  routed  to  charge  amplifiers,  which  converted 
the  charge  output  of  the  transducers  into  a  voltage,,  and 
then  to  the  wave  analyzers  for  display  as  RMS  values  and 
to  the  storage  oscilloscope,  where  the  waveforms  could  be 
compared . 


The  mass  flow 

rate 

through 

the 

test 

ine  was 

determined  using  a 

rotometer  that 

had 

been  calibrated 

against  a  gasometer, 

also 

known  as 

a  bell 

prover 

.  Since 

the  dynamic  response 

of 

the  line 

was 

altered 

by  its 

presence,  the  rotometer  was  removed  prior  to  each  test 
run.  The  rotometer  losses  were  small  enough  that  its 
removal  did  not  significantly  change  the  measured  mass 
flow  rate.  Variations  in  the  air  supply  during  the  test 
were  monitored  by  the  static  pressure  transducer. 

For  the  blocked  line  cases,  a  change  in  the  test  set 
up  was  required.  Since  the  driver  required  through  flow 
for  proper  operation,  a  bypass  valve  located  between  the 
driver  and  line  was  opened  effectively  turning  the  test 
line  into  a  blocked  branch  to  a  line  with  mean  flow. 

A  complete  list  of  the  major  instrumentation  used  can 


be  found  in  Appendix  D. 


Test  Designation 

A  5-digit  code  was  used  to  describe  each  test  case. 
The  first  digit  gave  the  line  number;  1  for  the  0.041,  2 
for  0.119  and  3  for  0.195  in.  inside  diameter  lines.  The 
second  digit  indicated  the  termination  employed;  0  for 
blocked,  1  and  2  for  0.0135  and  0.016  in  diameter  orifices 
and  3  for  open  lines.  The  last  three  digits  gave  the 
nominal  Reynolds  number  in  hundreds.  For  example,  the 
12020  case  was  the  0.041  in  test  line  fitted  with  a  0.0135 
in  orifice  termination  tested  at  a  Reynolds  number  of 
2000. 

Procedure 

Since  the  usable  frequency  range  of  the  driver 
assembly  was  limited  to  0  to  1000  Hz,  the  remaining  means 
to  vary  Q  involved  line  radius  and  kinematic  viscosity. 
For  a  given  line,  assuming  constant  absolute  viscosity  and 
perfect  gas  behavior,  the  range  was  dictated  by  line 
pressure  and  temperature,  taken  to  be  Pp  and  Ta 
respectively.  As  seen  in  Fig  3,  increasing  the  line 
pressure  in  a  given  line  will  allow  higher  £? values.  For 
turbulent  mean  flow  cases,  however,  increases  in  Pp 
yield  increases  in  attenuation  below  the  break  frequency, 
the  start  of  laminar  behavior.  This  point,  termed  Q g, 
is  then  shifted  to  correspondingly  higher  values  which  may 
exceed  the  maximum  possible  with  the  given  line.  Since  it 
is  of  particular  interest  to  investigate  the  response  of 
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.turbulent  flow  cases  in  the  mid  frequency  range,  orifice 
terminations  are  used  to  increase  the  PT  and  change  the 
Reynolds  number.  Reynolds  numbers  are  limited  by  a  50 
psig  line  pressure  constraint  imposed  by  the  static 
transducer  and  test  lines.  The  air  supply  was  manually 
regulated  to  provide  the  mass  flow  required  to  obtain  the 
desired  Reynolds  number.  The  rotometer  was  then  removed 
from  the  line  and  fluctuations  in  the  line  pressure  were 
monitored  by  the  static  pressure  transducer. 

The  frequency  of  the  sending  wave  analyzer,  and 
therefore  the  signal,  was  varied  from  20  to  1000  Hz  in  20 
Hz  increments.  At  each  increment,  the  frequency  of  the 
receiving  wave  analyzer  was  adjusted  to  provide  the 
maximum  voltage  reading.  The  delay  time  between  the 
sending  and  receiving  waveforms  was  read  off  the 
oscilloscope.  The  RMS  voltages  displayed  on  the  wave 
analyzers  and  delay  time  were  recorded  before  moving  to 
the  next  frequency.  These  data  became  part  of  the  input 
file  for  the  computer  program  and  were  used  to  generate 
the  experimental  gain  and  phase  shift. 


IV.  Results  and  Discussion 


Test  Conditions 

Table  III  summarizes  the  line  conf igurations  and  test 
conditions  that  were  used  in  an  attempt  to  examine  various 
regions  of  Fig  3.  The  range  of  nondimensional  frequency 
spanned  by  each  turbulent  test  case  is  also  presented. 
Referring  to  Fig  3,  it  is  seen  that  the  majority  of  the 
data  collected  is  within  the  mid  frequency  range. 

Blocked  Lines 

Fig  8-13  give  the  predicted  and  experimentally 
obtained  gain  and  phase  shift  plotted  versus  frequency  for 
the  three  blocked  line  cases.  These  figures,  along  with 
those  for  laminar  and  turbulent  mean  flow  cases,  are  found 
in  Appendix  A.  The  prediction  capaoilities  of  equations 
(6)  and  (7)  have  been  demonstrated  by  many  studies  to  be 
excellent  for  blocked  lines;  hence  providing  a  very 
reliable  means  to  assess  the  experimental  data  collection 
techniques  and  evaluate  the  line  model.  Fig  10-13  show 
good  agreement  between  theoretical  and  experimental  values 
for  the  two  larger  diameter  lines.  Fig  8  and  9,  however, 
show  greater  deviation  for  line  1,  the  0.041  in.  diameter 
line,  which  is  probably  attributable  to  the  increased 
difficulty  of  accurately  modeling  the  smaller  line 


dimensions . 


Table  III 
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Case 
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Table  III 


Laminar  Mean  Flow 


Fig  14-21  present  predicted  and  measured  gain  and 
phase  shift  for  the  cases  of  laminar  mean  flow  in  which 
the  flow  is  assumed  to  be  choked  in  the  orifice. 
Schlieren  techniques  were  applied  to  several  cases  to 
confirm  sonic  flow  at  the  the  orifice  exit. 

As  for  the  blocked  lines,  equations  (6)  and  (7)  were 
used  to  predict  response.  From  the  blocked  line  analysis, 
the  line  models  and  experiment  are  seen  to  be  reasonably 
accurate.  Therefore,  the  effect  of  the  terminating 
impedance  could  be  determined. 

For  comparison,  the  terminal  impedance  was  calculated 
using  two  different  models.  The  classic  expression  for 
the  steady  state  impedance  of  an  orifice  is  found  from 
Bernoulli's  equation  to  be 

ZL  =  [p  4P/  (2AQ2cd2)  ]0.5  +  j  o  (51) 

where  A0  is  the  orifice  area  and  C<j  is  the  discharge 
coefficient.  The  AC  value  is  shown  by  equation  (38)  to  be 
twice  the  valu^  obtained  in  equation  (51).  Fig  14-21 
display  response  using  terminal  impedances  obtained 
through  the  appropriate  reflection  coefficient,  equations 
(42)  and  (43)  and  the  AC  values  of  equation  (51). 

Althouch  the  AC  values  vary  from  approximately  equal 
to  more  than  twice  those  found  using  reflection 
coefficients,  the  response  for  the  two  impedances  is  seen 
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to  change  very  little.  This  is  due  to  the  ratio  of 
orifice  to  line  area.  As  this  ratio  goes  to  zero,  the 
line  will  appear  blocked.  For  lines  that  are  blocked,  or 
almost  blocked,  a  certain  impedance  value  exists  above 
which  response  changes  are  negligible. 

For  cases  in  which  the  pressure  ratio  across  the 
orifice  is  insufficient  to  choke  to  flow,  equation  (42) 
can  again  be  used  to  find  the  terminating  impedance  if  the 
reflection  coefficient  from  equation  (44)  is  used.  The 
impedances  for  these  cases  are  larger  than  those  for  the 
choked  cases  and  are  therefore  seen  to  be  in  better 
agreement  with  the  AC  impedance  values. 

Fig  22-33  present  experimental  gain  and  phase  shift, 
as  well  as  theoretical  predictions  using  the  terminal 
impedance  equal  to  the  AC  values  of  equation  (38)  and 
impedance  found  using  reflection  coefficients  in  equation 
(42),  for  the  unchoked  orifice  cases.  The  theory,  using 
either  terminal  impedance,  is  seen  in  Fig  22-29  to  predict 
higher  attenuation  than  that  observed  experimentally  for 
line  1.  The  theoretical  predictions  for  lines  2  and  3, 
diameters  of  0.119  and  0.195  in.  respectively,  are  seen  in 
Fig  30-33  to  be  in  much  better  agreement  with 
experimental  data.  The  overprediction  of  attenuation  for 
the  line  1  case  is  similar  to  that  seen  in  Fig  8-9  for  the 
same  line  in  the  blocked  configuration.  Since  only  the 
terminal  impedance  was  changed  in  going  from  the  blocked 


to  orifice  case,  it  appears  the  deviations  in  Fig  22-29 
are  related  to  the  problem  of  properly  modeling  line  1. 

Cases  of  laminar  flow  through  open  lines  offer  an 
opportunity  to  assess  the  effect  of  using  the  radiation 
impedance  as  given  by  equation  (41)  versus  the  value  given 
by  equation  (38)  which  becomes  identically  zero  for  open 
lines.  The  radiation  impedance  is  seen  to  be  quite  small 
at  all  dimensional  frequencies  of  this  study  and  as  shown 
by  Fig  34-37  may  be  neglected.  The  radiation  impedance 
will,  however,  be  used,  along  with  that  found  using 
reflection  coefficients  for  orifices,  for  the  remainder  of 
this  study. 

Turbulent  Mean  Flow 

The  blocked  and  laminar  cases  were  used  to  verify  the 
experiment,  line  model  and  terminating  impedance.  As  an 
initial  step  in  determining  turbulent  flow  response,  the 
components  of  the  constant  LRC  model  are  varied 
individually. 

Fig  39  shows  the  effect  of  varying  the  resistance  on 
the  gain  predictions  for  an  open  line.  AC  and  DC 
resistances,  from  equations  (23)  and  (31)  respectively, 
were  used  with  the  adiabatic  inertance  and  isothermal 
capacitance.  As  seen  by  Briski  (Ref  5),  use  of  the  AC 
resistance  yields  predicted  peak  gain  magnitudes  near 
those  seen  experimentally  but  occurring  at  lower 
frequencies.  This  shift  would  describe  the  response  of  a 


similiar  line  of  longer  length.  Since  the  response  of 
the  line  is  actually  dictated  by  the  product  of  the 
propagation  operator  and  line  length,  as  seen  in  equation 
(17),  an  increase  in  either  increases  the  effective  line 
length.  The  increased  attenuation  resulting  from  the  AC 
resistance  is  seen  to  delay  the  transition  to  the  high 
frequency  laminar  behavior  from  ,  in  this  case, 
approximately  340  Hz  for  the  DC  resistance  to  a  frequency 
beyond  the  1000  Hz  test  limit. 

Fig  39  illustrates  the  effect  of  using  either  the 
isothermal  or  the  turbulent  capacitance  with  the  AC 
resistance  and  adiabatic  inertance.  The  use  of  the 
smaller  turbulent  capacitance,  given  by  equation  (36), 
results  in  a  smaller  propagation  operator  which 
effectively  shortens  the  line.  The  peak  gain  magnitudes 
are  seen  to  be  reduced  very  little  and  now  occur  at 
slightly  higher  frequencies  than  the  experimental  data. 
Attenuation  remains  sufficiently  high  enough  to  inhibit 
transition  below  1000  Hz. 

Fig  40  shows  the  effect  of  varying  the  inertance. 
The  adiabatic  and  turbulent  inertance,  equations  (27)  and 
(32)  respectively,  were  used  with  the  AC  resistance  and 
turbulent  capacitance  to  obtain  the  theoretical  gain 
predictions.  For  this  case  of  Reynolds  number  equal  to 
10000,  the  turbulent  inertance  is  only  approximate ly  3  % 
greater  than  the  adiabatic  value  and  therefore  results  in 


only  a  slight  increase  in  the  effective  line  length.  The 
attenuation  shows  no  appreciable  change. 

The  combination  of  AC  resistance,  turbulent 
capacitance  and  turbulent  inertance  will  be  referred  to  as 
the  modified  constant  LRC  model.  The  modified  constant 
LRC  incorporates  the  imaginary  part  of  the  conductance, 
within  the  turbulent  capacitance,  but  neglects  the  real 
part.  Adding  the  real  part  of  equation  (35)  to  the 
modified  constant  LRC  model  yields  the  constant  LRCG 
mode  1 . 

Fig  55  shows  the  gain  predictions  for  the  modified 
constant  LRC  and  constant  LRCG  models  for  the  case 
presented  in  Fig  38-40.  The  constant  LRCG  model  predicts 
higher  attenuation  than  the  modified  constant  LRC  model 
and  a  slightly  shorter  effective  line  length. 

The  modified  constant  LRC,  constant  LRCG  and  laminar, 
from  equations  (6)  and  (7),  models  are  used  to  generate 
theoretical  gain  and  phase  shift  for  turbulent  mean  flow 
in  orifice  terminated  lines.  Fig  41-48,  and  open  lines. 
Fig  49-64. 

The  predictions  for  cases  involving  line  1  are  again 
seen  to  be  in  poor  agreement  with  experimental  data.  Fig 
41  shows  the  laminar  model  to  predict  a  higher  attenuation 
than  seen  experimentally.  This  is  similar  to  the  behavior 
of  the  laminar  cases  for  line  1  as  shown  in  Fig  22-29. 
The  low  attenuation  predicted  by  the  laminar  model  at  low 


frequencies  in  Fig  43-48  is  as  expected  (Ref  3-6)  but  is 
accompanied  by  a  most  unexpected  increase  in  the 
propagation  operator. 

This  behavior  may  be  due  to  several  factors;  the 
errors  in  modeling  line  1,  errors  in  determining  the 
terminal  impedance  and/or  the  possibility  of  appreciable 
roughness  in  the  line.  The  friction  factor  given  by 
equation  (22)  is  applicable  to  turbulent  flow  in  smooth 
lines.  If  the  lines  are  not  smooth,  the  actual  friction 
factor  is  larger  with  the  deviation  from  the  smooth  values 
increasing  with  Reynolds  number.  Since  the  computer 
program  assumes  smooth  lines  to  calculate  pressure  drops 
along  the  line,  the  theory  would  not  yield  accurate 
predictions  if  the  line  were  sufficiently  rough. 

Due  to  these  problems,  line  1  can  not  be  used  to 
assess,  with  adequate  confidence,  the  prediction 
capabilities  of  the  constant  LRCG  and  modified  constant 
LRC  models. 

Fig  49-64  show  experimental  theoretical  gain  and  phase 
shift  for  open  line  cases  for  lines  2  and  3,  diameters  of 
0.119  and  0.195  in.  The  laminar  model  predictions  are  as 
expected  indicating  little  line  model  error  and 
insignificant  roughness  effects.  Since  the  roughness 
effects  are  a  function  of  the  relative  roughness,  defined 
as  the  line  roughness  divided  by  the  line  diameter,  the 
larger  lines  would  exhibit  less  increase  in  friction 
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factor  than  the  smaller  line  with  comparable  line 
roughness . 

Fig  49-64  show  the  modified  constant  LRC  model  to 
predict  less  attenuation  than  the  constant  LRCG  model. 
The  laminar  model  is  shown  to  correctly  predict  the 
effective  line  length  while  the  modified  constant  LRC  and 
constant  LRCG  models  slightly  underpredict  the  propagation 
operator . 

Examination  of  the  assumptions  and  simplifications 
made  in  the  formulation  of  both  the  modified  constant  LRC 
and  constant  LRCG  models  reveals  some  possible 
explanations  for  the  degraded  prediction  capabilities  for 
certain  cases.  To  obtain  equation  (35)  from  equation 
(34),  the  term  (R^/of  )  was  assumed  to  be  much  smaller 
than  the  turbulent  inertance.  However,  in  cases  of  large 
Reynolds  number  and/or  small  line  diameter,  the  AC 
resistance  can  become  large  enough  to  violate  this 
assumption.  Also,  this  assumption  does  not  hold  for  very 
low  frequencies.  The  inclusion  of  this  term  will  tend  to 
lower  the  conductance  resulting  in  lower  attenuation. 
Comparison  of  the  two  larger  lines  tested  at  a  Reynolds 
number  of  5000,  Fig  51  and  59,  indicate  increasing 
overprediction  of  the  attenuation  as  the  line  size  is 
decreased.  Fig  49-58  are  seen  to  indicate  similar 
increases  in  attenuation  error  due  to  increasing  Reynolds 
number.  The  0.195  in  diameter  line  is  apparently  large 


enough  that  the  term  in  question  may  be  ignored  at  the 
frequencies  and  Reynolds  numbers  tested.  Incorporation  of 
this  term  should  enhance  the  prediction  capability  of  both 
models  but  stray  from  the  original  intent  of  keeping  the 
components  frequency  independent. 
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V.  Conclusions 


1.  The  impedance  of  orifices  with  and  without  choked  flow 
could  be  adequately  determined  using  reflection 
coefficients.  The  small  radiation  impedance  associated 
with  the  open  lines  tested  could  be  neglected  for 
frequencies  less  than  1000  Hz. 

2.  The  laminar  model,  using  the  Krishnaiyer  and  Lechner 
equations,  gave  accurate  predictions  for  blocked  test 
lines  and  lines  tested  with  laminar  mean  flow.  Attempts 
to  apply  this  model  to  turbulent  mean  flow  cases  resulted 
in  good  agreement  with  experimental  data  at  higher 
frequencies,  but  significant  underprediction  of  the 
attenuation  at  lower  frequencies. 

3.  The  modified  constant  LRC  model,  consisting  of  the 
turbulent  inertance,  AC  resistance  and  turbulent 
capacitance,  was  seen  to  predict  higher  attenuation  than 

>  the  laminar  model  for  turbulent  mean  flow  cases  at  lower 

frequencies.  As  a  result,  deviation  between  theoretical 
and  experimental  peak  gain  magnitudes  was  reduced,  for  the 
cases  tested,  by  at  least  50  %  while  introducing  less  than 
3  %  error  in  the  frequency  at  which  the  peaks  occur. 

4.  The  constant  LRCG  model,  which  adds  the  real  part  of 
the  conductance  to  the  modified  constant  LRCG  model. 
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further  increases  the  attenuation  resulting  in  peak  gain 
magnitudes  approximately  3  db  less  than  those  predicted  by 
the  modified  constant  LRC  model  and,  in  most  cases,  well 
below  the  experimental  values. 

5)  The  modified  constant  LRC  and  constant  LRCG  models 
account  for  the  effects  of  Reynolds  number,,  but  not 
frequency,  on  the  individual  components.  The  applicable 
frequency  range  of  both  models  is  limited  by  the 
combination  of  Reynolds  number,  density  and  line  diameter 
due  to  certain  assumptions  made  to  eliminate  frequency 


dependence . 


VI.  Recommendations 


1.  The  effect  of  including  the  frequency  dependent 
terms  of  the  conductance  should  be  investigated.  The 
improvement  in  prediction  capability  should  be  weighed 
against  the  increased  complexity  introduced. 

2.  Experimental  data  points  should  be  extended  into 
the  low  frequency  range.  Present  instrumentation,  in 
particular  the  wave  analyzers,  prevents  accurate  readings 
at  the  extremely  low  frequencies  required  to  explore  this 
region  with  the  lines  used  this  study.  Therefore,  using 
current  instrumentation,  smaller  line  diameters  must  be 
considered  to  lower  the  nondimensional  frequency. 

3.  Automated  data  collection  equipment  should  be 
used  to  increase  accuracy,  facilitate  low  frequency 
measurements,  standardize  experimental  procedures  and 
significantly  reduce  time  requirements  while  providing 
more  data  points.  An  experimental  data  base  could  then  be 
established,  using  reliable  and  repeatable  results 
obtained  through  standardized  test  procedures ,  that  future 
studies  could  draw  upon  and  contribute  to. 

4.  The  effect  of  line  roughness  on  the  response  of 
turbulent  flow  lines  should  be  investigated  in  greater 
detail . 
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Appendix  A 


□  a 


Figure  9  Experimental  and  Theoretical  Phase  vs  Frequency  for  Case  10000 


imental  and  Theoretical  Gain  vs  Frequency  for  Case  20000 


Figure  12  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  30000 
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Theory  with  Z  from  Eq  (38) 


Figure  18  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  31010 


Theory  with  Z  from  Eq  (38) 


Theory  with  Z  from  Eq  (38) 


Theory  with  Z  from  Eq  (39) 
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Figure  22  Experimental  and  Theoretical  Gain  vs  Frequency  f 


Theory  with  Z  from  Eq  (38) 


Figure  23  Experimental  and  Theoretical  Phase  vs  Frequency  for  Case  11010 


Laminar  Theory  with  ZT  from  Eq  (39) 


Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  11020 


Theory  with  Z  from  Eq  (39) 


Figure  25  Experimental  and  Theoretical  Phase  vs  Frequency  for  Case  11020 


Theory  with  Z  from  Eq  (39) 


Figure  26  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  12010 


Theory  with  Z  from  Eq  (39) 


Figure  27  Experimental  and  Theoretical  Phase  vs  Frequency  for  Case  12010 
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ental  and  Theoretical  Gain  vs  Frequency  for  Case  12020 


Theory  with  Z  from  Eq  (38) 
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Figure  33  Experimental  and  Theoretical  Phase  vs  Frequency  for  Case  31005 


mental  and  Theoretical  Gain  vs  Frequency  for  Case  23020 


imental  and  Theoretical  Gain  vs  Frequency  for  Case  33020 
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Figure  37  Experimental  and  Theoretical  Phase  vs  Frequency  for  Case  33020 


Figure  38  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  23100 


igure  39  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  23100 


Figure  41  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  11030 
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Figure  45  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  12050 
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Figure  51  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  23050 
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Figure  55  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  23100 
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Appendix  B 


Orifice  Models 

A  number  of  different  methods  were  considered  to 
determine  the  impedance  of  an  orifice  with  through  flow. 
The  gain  predictions  of  the  various  models  are  shown  for 
case  22020  in  Fig  65. 

Model  1 

Since  the  orifice  diameter  d  is  much  smaller  than  the 
line  diameter  D,  blocked  line  impedance  was  assumed 

ZL  =  ob  +  j  0  (52) 

Model  2 

Equation  (51)  was  used  assuming  constant  density  equal 
to  that  just  before  the  orifice,  negligible  length  and 
discharge  coefficient  of  0.6. 

Model  3 

The  AC  impedance  of  equation  (51),  seen  to  be  twice 
the  value  obtained  from  equation  (51),  is  used  with  the 
same  assumptions  as  in  Model  2. 

Model  4 

Model  2  was  modified  to  use  mean  density  by  assuming 
the  pressure  in  the  orifice  to  be  half  the  gauge  pressure 
before  the  orifice. 


Model  5 


The  same  average  density  was  used  in  conjunction  with 
the  AC  impedance  resulting  in  twice  the  impedance  of  Model 

4. 

Model  6 

The  steady  state  DC  impedance  of  Model  2  is  modified 
by  assuming  a  sonic  density  in  the  orifice.  The  sonic 
density  is  found  through  the  use  of  continuity  to  be 

sonic  =  M/ ( A0 {  VkRgTa ) )  (53) 

where  the  Mach  number  has  been  taken  as  unity.  Using  the 
sonic  density,  the  DC  impedance  can  be  expressed  as 


ZL  =  pl/Ql 


(2 (0.6) 2A  3(kR  t a)0.5 


(54) 


Model  7 

The  AC  impedance  was  calculated  using  as  in  Model  6,  a 
Mach  number  of  1  in  the  orifice  as  in  Model  6. 


Model  0 

""  • 

This  model  assumes  the  flow  to  be  isentropic  and 
choked  with  the  A*  value  taken  to  be  equal  to  the 
orifice  area.  Since  d<<D,  the  flow  in  the  line  is  very 
small  and  the  properties  upstream  of  the  orifice  are  taken 
to  be  the  stagnation  values.  With  these  assumptions,  the 
DC  impedance  becomes 

ZL  P  PL 
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Model  9 


With  the  same  assumptions  as  in  Model  8,  the  AC 
impedance  is  again  twice  the  DC  value. 

Model  10 

Model  8  is  used  with  the  density  take  to  be  the 
average  value  previously  mentioned. 

Model  11 

The  AC  impedance  of  Model  9  is  calculated  using  the 
average  density  value. 

Model  12 

Model  8  is  used  with  the  sonic  density  to  obtain  a  DC 
impedance . 

Model  13 

The  sonic  density  is  substituted  for  the  constant 
density  in  Model  9. 

Model  14 

Funk,  Wood  and  Chao  (Ref  18)  present  an  extension  to 
the  DC  impedance  of  Model  2  by  considering  the  axial 
length.  The  length  of  the  orifice  is  taken  to  be  0.06 
which  is  approximately  4  diameters  in  this  case.  The  flow 
is  assumed  to  be  incompressible  with  a  discharge 
coefficient  of  0.6.  The  friction  factor  is  found  using  a 
Reynolds  number  of  19000  in  the  orifice.  The  resistance 


contributed  by  the  orifice  axial  length  is  seen  to 
increase  the  impedance  by  less  than  5  percent. 


Model  15 

The  above  method  is  applied  with  a  correction  to  the 
discharge  coefficient  for  compressible  flow.  Shapiro  (Ref 
19)  shows  the  discharge  coefficient  for  sharp-edged 
orifices  to  vary  in  a  non-linear  fashion  with  the  ratio  of 
the  downstream  to  upstream  pressures.  The  discharge 
coefficient  ranges  from  0.6,  as  is  assumed  for 

incompressible  flow,  to  approximately  0.85  as  the  pressure 
ratio  goes  from  1.0  to  0.2.  Greater  pressure 

differentials  yield  negligible  increases  beyond  the  0.85 
value.  For  this  test  case,  the  pressure  ratio  is  0.28 
giving  a  discharge  coefficient  of  0.82.  Model  14  is  then 
calculated  using  this  value. 

On  close  examination  under  a  microscope,  the  orifices 
used  for  this  study  were  seen  to  be  far  different  from 

sharp-edged  or  even  square-edged  models.  Since  the 

problem  of  compressible  flow  in  an  orifice  is  very 
difficult  even  when  dealing  with  ideal  models,  alternative 
approaches  to  determining  the  impedance  that  do  not 

require  the  orifice  density  and  discharge  coefficient  are 
sought . 

Model  16 

The  use  of  reflection  coefficients  is  outlined  in 
Section  II.  When  a  pressure  wave  propagating  along  a  line 


[»] 


encounters  an  area  change,  as  shown  in  Fig  4,  a  reflected 
wave  is  generated,  the  size  and  phase  of  which  is 
determined  by  the  area  change.  In  the  case  of  a  blocked 
line,  the  original  and  reflected  waves  are  identical.  The 
reflection  coefficient  for  this  case  is  then  defined  to  be 
1.  In  the  other  extreme  of  an  open  line,  the  downstream 
area  is  infinite  and  the  reflected  pressure  wave  is  equal 
and  180  degrees  out  of  phase  with  the  original.  The 
corresponding  reflection  coefficient  is  -1.  For  the  case 
where  there  is  no  area  change,  no  reflected  waves  are 
generated  and  the  reflection  coefficient  becomes  0.  This 
is  also  seen  to  be  the  case  when  sonic  flow  exists  as  the 
reflected  waves  cannot  propagate  up  stream  against  the 
mean  flow. 

To  find  the  reflection  coefficients  for  orifices, the 
cross-sectional  area  at  the  end  of  the  line  as  seen  by  the 
pressure  wave  can  be  considered  as  three  distinct 
regions.  The  orifice  is  split  into  two  areas,  a  core  area 
in  which  the  flow  is  sonic,  given  a  sufficient  pressure 
ratio,  and  the  remaining  annular  area  which  is  open  to  the 
ambient  conditions.  The  annular  area  around  the  orifice 
acts  as  a  blocked  line.  The  reflection  coefficient  of 
each  region  is  then  multiplied  by  its  area  to  determine 
its  contribution  to  an  average  reflection  coefficient 
which  can  be  used  in  equation  (42)  to  find  the 


For  unchoked  cases,  no  core  area  exists  and  the  model 
is  reduced  to  two  regions.  For  choked  cases,  it  was 
believed  that  the  vena-contracta  area  would  be  the  core 
area.  An  attempt  to  determine  this  area  using  Schlieren 
flow  visualization  indicated  a  region  of  supersonic  flow 
at  the  orifice  exit.  The  flow  pattern  suggests  the  sonic 
flow  had  fully  filled  the  orifice  area  before  the  exit 
plane  thus  allowing  the  use  of  a  two  region  model  where 
the  core  area  is  the  entire  orifice  are. 

Model  16 

As  the  orifice  diameter  approaches  the  line  diameter, 
the  impedance  tends  towards  that  for  an  open  line.  The 
limiting  case  of  the  impedance  equal  to  zero  is  shown  for 
completeness . 
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APPENDIX  C 


Computer  Program 


PROGRAM  MSB 

DIMENSION  DBT ( 300 ) ,RVT(100) ,DI (100) , AD (100) ,AR(100) , 
1CV(100) , OMT (100) , CNA (100) ,GMA(100) ,AGM(100) ,FN(100) , 
2RN ( 100 ) , RN (100) ,GN(100) ,ALN(100) ,CN(100) ,AN(100) , 

3BTN ( 100 ) , AMC (100) ,AZRN(100) ,BZRN(100) ,AZIN(100) , 

4BZIN ( 100) , DC ( 100 ) ,DG(100) ,DD(100) ,P(100) ,RHO(100) , 
5ANU (100) , CA ( 100 ) ,QTCA(100) ,OMG(300) ,RTP(100) ,GP(100) , 
6BETA ( 100 ) , AZOT ( 100 ) ,BZOT(100) ,REY(100) ,Q(100) , 

7BET (100) ,TURB (100) 

DIMENSION  DB ( 300 ) ,OMGX(300) ,OMGP(300) ,PHASE(300) , 

1PH (300) ,CASE ( 100 ) 

DATA  PI/3.1415926/ 

DATA  TPI/6. 2831853/ 

C 

OPEN ( 1 1 , FILE= 1 D10 00  1  ) 

REWIND  11 

OPEN ( 12 , F ILE= 1  AC ROUT ' ) 

REWIND  12 

CALL  PLOTS (0. ,0. ,9) 

C  TAPE  11  IS  THE  INPUT  FILE  CONTAINING  THE  FOLLOWING 
C  ICAS  IS  THE  CASE  NUMBER 

C  OP  INDICATES  OPEN  LINE;  1.0  FOR  OPEN  LINE,  ELSE  0.0 

C  TF  IS  AMBIENT  TEMPERATURE  IN  DEGREES  F 

C  PG  IS  THE  BAROMETRIC  PRESSURE  (PSIA) 

C  AMU  IS  THE  DYNAMIC  VISCOSITY  (PSI-SEC) 

C  RE  IS  THE  GAS  CONSTANT  FOR  AIR  (IN  /SEC  /  R) 

C  GAM  IS  RATIO  OF  SPECIFIC  HEATS 

C  SIG  IS  SQUARE  ROOT  OF  PRANDTL  NUMBER 

C  N  IS  NUMBER  OF  LINE  SEGMENTS 

C  CT  IS  THE  TURBULENT  CAPAC ITANCE/ ADIABATIC  CAPACITANCE 

C  TL  IS  THE  INERTANCE  CONSTANT 

C  FL  IS  THE  MASS  FLOW  RATE  (LBM/SEC) 

C  D1  IS  OUTSIDE  DIAMETER  OF  LINE  (IN) 

C  PT  IS  STATIC  LINE  PRESSURE  (PSIG) 

C  ZEND  IS  REAL  PART  OF  TERMINAL  IMPEDANCE  (PSI/CIS) 

C  ZENDI  IS  IMAGINARY  TERMINAL  IMPEDANCE  (PSI/CIS) 

C  DI  IS  THE  LINE  SEGMENT  DIAMETER  (IN) 

C  AD  IS  THE  LINE  SEGMENT  LENGTH  (IN) 

C  THE  FOLLOWING  ARE  USED  BY  THE  CALCOMP  PLOTTER 

C  NPTS  IS  NUMBER  OF  EXPERIMENTAL  POINTS  INPUT 

C  LSMB  DESIGNATES  PLOTTER  SYMBOL  (SEE  USERS  GUIDE) 

C  FREQ  IS  EXPERIMENTAL  FREQUENCY  (HZ) 

C  PS  IS  SENDING  PRESSURE,  EXPERIMENTAL 

C  PR  IS  RECEIVING  PRESSURE,  EXPERIMENTAL 

C  PHD  IS  EXPERIMENTAL  PHASE  DELAY  IN  MSEC 

110  READ (11,*) ICAS, OP 
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READ ( 1 1 , * ) TF , PG , AMU , RE , GAM , S IG 

READ ( 1 1 ,  * )  N 

READ (11,*) CT , TL 

READ (11,*) FL 

READ (11,*) D1,PT 

READ (11,*) ZEND , ZENDI 

DO  898  1=1, N 

READ ( 11 , * ) DI ( I ) , AD ( I ) 

CONTINUE 
NPG  =  0 
ICT  =  0 
IND  =  0 
NNN=0 
M=0 

IF (OP . EQ . 1 • )  GOTO  900 

REYT=1.22999E+06*FL/D1 

CONTINUE 

PT=PT+PG 

I F ( REYT . LT .2300.)  THEN 

DP=PT-SQRT (PT**2 .+2 . 27E-10*REYT* (TF+460.)/(Dl**3.)) 
ELSE 

DP=PT-SQRT (PT**2 .+1 . 12226E-12*REYT**1.75* (TF+460. ) / 
1 (Dl**3 . ) ) 

END  IF 

P (N) =PT+DP-PG 
PF=P (N) 

DO  899  1=1, N-2 

IF (OP .EQ . 1 . )  GOTO  895 

REYT=1 . 22999E+06*FL/DI (10) 

CONTINUE 

PF=PF+PG 

K=N+1-I 

IF (REYT. LT. 2300 . )  THEN 

DP=PF-SQRT (PF**2.+1 . 13 5E-10*AD ( K ) *REYT* (TF+460. ) / 

1 (Dl (K) **3  . )  ) 

ELSE 

DP  =  PF-SQRT (PF**2.  +  5.6 1 12  8E-1 3  *REYT*  *  1 . 75*AD ( K ) * 

1 (TF  +  460 . ) / (D1(K)**3.)  ) 

END  IF 

P(N-I) =PF+DP-PG 
PF=P (N-I ) 

PEND=P (2) 

I F (OP . EQ . 1 .  )  PEND=0 . 0 
CONTINUE 
P ( 2 ) =P ( 5 ) 

P(l) =P(2) 

GOTO  903 
CONTINUE 
TOTL=0 .0 
DO  901  1=3, N 
TOTL=TOTL+AD( I) 

CONTINUE 
TOTL  =  TOTL+  2 . 
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REYT=PT* (PT/2.+PG) * (DI (N) *  *  3 . ) / 32 . / AMU/ AMU 
REYT=REYT/TOTL/RE/ (TF  +  460  . ) 

I F ( RE YT . LT .2300.)  GOTO  902 
REYT = ( REYT* 32./.1582) **{4./7.) 

2  CONTINUE 
GOTO  897 

3  CONTINUE 
TOTL=0 . 0 
DO  2  1=3, N 
CONTINUE 

RHO(N) = { PG+P (N) ) / (RE* (TF+460. ) ) 

REYT=  (P  (N)  -PEND)  *  (  (P  (N)  +PEND)  /  2  . +PG)  *  (DI  (N)  * *3  . )  /  32  .  / 
1 AMU / AMU / TOTL/ RE / (TF  +  460. ) 

I F ( REYT . LT .2300. )  GOTO  3 
REYT= ( REYT*  3  2./.1582)**(4./7.) 

QM=REYT*PI*DI (N) *AMU/4. 

DO  23  1=1, N 
PBR=P ( I ) +PG 
TBR=TF  +  460 . 

RHO(I) =PBR/ ( RE*TBR) 

ANU(I) =AMU/RHO ( I ) 

CA ( I ) =SQRT ( PBR*GAM/ RHO ( I ) ) 

AR (I) =PI*DI (I) *DI (I) /4. 

CV(I) = (8.*PI*ANU(I) ) / AR ( I ) 

OMT(I) =CV(I) /SIG/SIG 

CNA ( I ) = ( 8 . *PI *AMU ) / AR ( I ) / AR ( I ) 

Q ( I ) =QM/RHO ( I ) 

CONTINUE 
Q(l) =0.0 
Q (2) -0.0 
GMN I = . 5  * (GAM-1. ) 

DO  26  1  =  1, N 

REY ( I )  =  ( 4 . *  RHO ( I ) *Q ( I ) /  (PI*DI ( I) *AMU) 

PBR=P ( I ) +PG 
TEMP=GMNI/ (GAM*PBR) 

GMA(I) =TEMP*AR(I) 

AGM(I) =AR(I) / (GAM*PBR) 

QTCA(I) =0.25*CA(I) 

RVT ( I ) = . 3 164* AMU/ 2 . / AR ( I ) / DI (I) /DI ( I ) * (REY ( I ) ** . 75) 
RVT(I) =  RVT ( I ) *1.75 
FN ( I ) =QTCA ( I ) / AD ( I ) 

NST=  1 
DW=0 . 

Y=1 . 

DO  80  J=1,NST 

M  =  M+ 1 

Y  =  Y  +  DW 

W=TPI*Y 

DO  27  1  =  1, N 

ARG=.5*SQRT(W/CV(I) ) 

RN (I) =CNA ( I ) * ( . 375+ARG+ ( .375/ (4.*ARG) ) ) 

DO  28  1  =  1, N 

DC (I) = . 25+SQRT (W/OMT ( I ) + . 1 25*SQRT (OMT ( I ) /W) 
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DG ( I ) =SQRT (W/ OMT ( I ) ) - . 12 5 *SQRT (OMT ( I ) /W) 

DD ( I ) =DC ( I ) *  DC ( I ) +DG ( I ) *DG ( I ) 

GN (I) =W* (GAM-1 . ) *AGM (I) *DG (I) /DD  (I) 

DO  29  1=1, N 
ARG= . 5*SQRT (CV ( I ) /W) 

ALN(I) =RHO(I) * (l.+ARG- (ARG* (15.*CV(I) / (W*64. ) ) ) ) / 
1AR(I) 

TEMP=GMNI/W 
DO  30  1  =  1, N 
TURB(I) =0 

CN ( I ) =AGM ( I )  *  ( 1 .  +  ( (GAM-1. ) *DC (I) /DD(I) ) ) 

TEMP=-W*W 
DO  31  1=1, N 

TEM1=RN ( I ) *GN ( I ) +TEMP*ALN ( 1 ) *CN ( I ) 

TEM2=W* ( RN ( I ) *CN ( I ) +GN(I) ALN(I) ) 

CALL  RTCMP ( ARG1 , ARG2 ,TEM1 ,TEM2 ) 

AN ( I ) = ARG1 
BTN ( I ) = ARG2 

C  NEXT  2  TERMS  ARE  DETERMINED  BY  THE  TURBULENT  MODEL 

C  TEM1  IS  THE  REAL  PART  OF  ZY 

C  TEM2  IS  IMAGINARY  PART  OF  ZY 

C  THE  TERMS  SHOWN  ARE  FOR  THE  CONSTANT  LRCG  MODEL 
C  FOR  MODIFIED  CONSTANT  LRC  USE 

TEM1=— W*W*TL*RHO ( I ) / AR ( I ) *CT* AGM ( I ) 

TEM2=W*RVT ( I ) *CT*AGM ( I ) 

************************ 

TEM1=RVT ( I ) *RVT ( I ) * (GAM-1 . ) *CT* AGM ( I ) *AR ( I ) /RHO  ( I ) / 
1TL/TL-W*W*TL*RH0 ( I ) / AR ( I ) *CT*AGM ( I ) 

TEM2  =W*RVT ( I ) *CT* AGM ( I ) +W*RVT ( I ) * (GAM-1 . ) *AGM ( I ) /TL 
CALL  RTCMP ( ARG1 , ARG2 , TEM1 , TEM2 ) 

TAN = ARG 1 
TBTN=ARG2 

I F ( RE Y ( I ) . LT .2300. )  GOTO  31 
I F ( AN ( I ) . GE . TAN )  GOTO  31 
AN ( I ) =TAN 
BTN ( I ) =TBTN 
TURB(I) =1 
AMC ( I ) = TP I /BTN ( I ) 

OX=AR ( 10 ) *W/ANU (10) /PI 
OY=AR ( 10 ) *CA ( 10 ) *AN ( 10 ) / PI / ANU ( 10 ) 

CALCULATE  Z0 
DO  32  1  =  1, N 
TEM1 =W*ALN (I) 

TEM2=W*CN( I) 

TEM3  =RN ( I ) 

TEM4=GN ( I ) 

CALL  CMPDV ( ARG1 , ARG 2 , TEM3 , TEM1 , TEM4  ,  TFM2 ) 

CALL  RTCMP ( AZRNI , BZRNI , ARG1 , ARG 2 ) 

AZRN(I) =AZRNI 
BZRN(I) =BZRNI 
IF(TURB(I) .EQ.0)  GOTO  32 
TEM1 =RHO ( I ) * ( P ( I ) +PG) / AR ( I ) / AR ( I ) 


TEM2=-RVT ( I ) * (P(I) +PG) /W/AR(I) 

CALL  RTCMP ( ARG1 , ARG2 , TEM1 , TEM2 ) 

AZRN (I) = ARG1 
BZRN(I) =ARG2 
CONTINUE 

CALCULATE  Z  IN  1 
1  =  1 

TEMP=AN ( I ) * AD ( I ) 

I F (TEMP . GT .88.)  GOTO  80 
ARGl=COSH (TEMP) 

ARG2=SINH (TEMP) 

TEM5=BTN ( I ) *AD(I) 

TEMP=COS (TEM5) 

TEM1=ARG1*TEMP 
TEM3  =ARG2*TEMP 
TEMP=SIN (TEM5 ) 

TEM2=ARG2*TEMP 

TEM4=ARG1*TEMP 

CALL  CMPDV (ARG1 ,  ARG2 ,TEM1 ,TEM2 ,TEM3 ,TEM4) 

TEM1=AZRN ( I ) 

TEM2=BZRN ( I ) 

CALL  CMPMP (TEM3 , TEM4 , TEM1 , TEM2 , ARG1 , ARG2 ) 

AZIN (I) =TEM3 
BZIN(I) =TEM4 
CALCULATE  Z  IN  3 
1  =  3 

TEM1=AZRN ( I ) 

TEM2=BZRN ( I ) 

IF (OP . EQ . 1 . )  THEN 
TEM3  =  0 . 0 

TEM4=-W*0 . 4*RHO ( I ) *DI ( I ) / AR ( I ) 

ELSE 

TEM3=ZEND 
TEM4=ZENDI 
END  IF 
TEM5  =AN ( I ) 

TEM6  = AD ( I ) 

TEM7=BTN ( I ) 

CALL  CALZ IN ( AARG , BARG , TEM1 , TEM2 , TEM3 , TEM4 , TEM5 , TEM6 
1TEM7 ) 

AZIN(I) = AARG 
BZ IN ( I ) =BARG 
CALCULATE  Z  IN  4  AND  5 
DO  34  1=4,5 
TEM1=AZRN( I) 

TEM2=BZRN ( I ) 

TEM3=AZ IN ( I — 1 ) 

TEM4=BZ IN ( 1-1 ) 

TEM5=AN ( I ) 

TEM6  =AD ( I ) 

TEM7=BTN ( I ) 

CALL  CALZ IN (AARG, BARG, TEM1 ,TEM2 ,TEM3 ,TEM4 ,TEM5,TEM6 
1TEM7 ) 


AZIN (I) =AARG 
34  BZ IN ( I ) =BARG 

C  CALCULATE  Z  IN  2 
1  =  2 

TEM1=AZRN ( I ) 

TEM2=BZRN ( I ) 

TEM3=AZIN ( I — 1 ) 

TEM4=BZ IN { I -1 ) 

TEM5=AN ( I ) 

TEM6=AD ( I ) 

TEM7  =BTN { I ) 

CALL  CALZ IN ( AARG , BARG , TEM1 , TEM2 , TEM3 , TEM4 , TEM5 , TEM6 , 
1TEM7) 

AZ IN ( I ) =AARG 
BZIN(I) =BARG 

C  CALCULATE  RECEIVING  Z  FOR  LINE  6 
1  =  6 

TEM1=AZIN ( I — 1 ) 

TEM2=BZ IN ( 1-1 ) 

TEM3=AZ IN ( 1-4 ) 

TEM4=BZ IN ( 1-4 ) 

CALL  ZEBRA ( AZOTI , BZOT I , TEM1 , TEM2 , TEM3 , TEM4 ) 

AZOT (I) =AZOTI 
BZOT (I) =BZOTI 
C  CALCULATE  Z  IN  6 

TEM1=AZRN ( I ) 

TEM2=BZRN ( I ) 

TEM3=AN ( I ) 

TEM4=AD(I) 

TEM5=BTN ( I ) 

CALL  CALZ IN (AZ INI ,BZ INI ,TEM1,TEM2, AZOTI ,BZOTI , TEM3 , 
1TEM4 , TEM5 ) 

AZIN(I) =AZ INI 
BZIN(I) =BZ INI 

C  CALCULATE  Z  IN  7,8,9,10 

DO  39  1=7,10 
TEM1 =AZRN ( I ) 

TEM2  =BZRN ( I ) 

TEM3  =AZ IN ( 1-1 ) 

TEM4  =  BZIN ( I  —  1 ) 

TEM5  =AN ( I ) 

TEM6  =AD ( I ) 

TEM7=BTN ( I ) 

CALL  CALZ IN (AZ INI , BZINI ,TEM1 ,TEM2 ,TEM3 ,TEM4 ,TEM5 , 
1TEM6 ,TEM7) 

AZIN(I) = AZ INI 
39  BZIN(I) =BZINI 
C  CALCULATE  P5/P6 
1=6 

TEMP=BTN ( I ) *AD ( I ) 

CSBI l=COS (TEMP) 

SNBI1=SIN (TEMP) 

TEMP=AN ( I ) *AD ( I ) 
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ARG1=C0SH (TEMP) 

ARG2=SINH (TEMP) 

TEM1=ARG1*CSBI1 
TEM2=ARG2*SNBI 1 
AZOTI =AZOT ( I ) 

BZOTI=BZOT ( I ) 

AZ INI = AZ IN ( I ) 

BZ INI =BZ IN ( I ) 

CALL  CMPDV ( TEM7 , TEM8 , AZOTI , BZOTI ,AZINI , BZ INI ) 
CALL  CMPMP (TEM3 , TEM4 , TEM7 , TEM8 , TEM1 , TEM2 ) 
TEM5=ARG2*CSBI 1 
TEM6=ARG1*SNBI1 
TEM9  =AZRN ( I ) 

TEM10=BZRN ( I ) 

CALL  CMPDV (TEM7 , TEM8 , AZOTI , BZOTI , TEM9 , TEM10 ) 

CALL  CMPMP (TEM1 , TEM2 , TEM7 , TEM8 , TEM5 , TEM6 ) 

TEM1=TEM3-TEM1 

TEM2=TEM4-TEM2 

TEMP=TEM1*TEM1 +TEM2*TEM2 

BETA (I) =ATAN2 (TEM2 , TEM1 ) 

BET (I) = (180. /PI) *BETA ( I ) 

RTP ( 1-1) =  SQRT (TEMP) 

GP(I-l) =20.*ALOG10 (RTP(I-l) ) 

CALCULATE  P6/P7,  P7/P8,  P8/P9,  AND  P9/P10 
DO  43  1=7,10 
TEMP=BTN ( I ) *AD ( I ) 

CSBIl=COS (TEMP) 

SNBI1=SIN (TEMP) 

TEMP= AN ( I ) * AD ( I ) 

ARGl=COSH (TEMP) 

ARG2=SINH (TEMP) 

TEM1=ARG1*CSBI1 

TEM2=ARG2*SNBI1 

AZINI1=AZIN(I-1) 

BZINI1=BZIN ( 1-1) 

AZINI =AZIN ( I ) 

BZ INI =BZ IN ( I ) 

CALL  CMPDV (TEM7 ,TEM8 , AZINI 1 ,BZ INI 1, AZINI ,BZ IN I) 
CALL  CMPMP ( TEM3 , TEM4 , TEM7 , TEM8 , TEM1 , TEM2 ) 
TEM5=ARG2*CSBI 1 
TEM6=ARG1*SNBI 1 
TEM9  = AZRN ( I ) 

TEM10  =BZRN ( I ) 

CALL  CMPDV (TEM7 , TEM8 , AZ INI  1 , BZ INI 1 , TEM9 , TEM10 ) 

CALL  CMPMP (TEM1 , TEM2 , TEM7 , TEM8 , TEM5 , TEM6 ) 

TEM1-TEM3-TEM1 

TEM2  =  TEM4  =TEM2 

TEMP=TEM1*TEM1 +TEM2*TEM2 

BETA (I) =ATAN2 (TEM2 ,TEM1 ) 

BET (I) = (180./PI) *BETA(I) 

RTP(I-l) =SQRT (TEMP) 

GP(I-l) =20.*ALOG10 (RTP(I-l) ) 

RTPT=RTP ( 6 ) *RTP ( 7 ) *RTP(8) *RTP(9) *RTP(5) 


GPT=20 . * ALOGIO { RTPT ) 

BETAR=BETA (6) +BETA ( 7 ) +BETA ( 8 ) +BETA ( 9 ) +BETA(10) 

35  IF (BETAR.LE.O. )  GOTO  36 
BETAR=BETAR-TPI 

GOTO  35 

36  IF (BETAR.GT.-TPI)  GOTO  37 
BETAR=BETAR+TPI 

GOTO  36 

37  BETAD= (180./PI) *BETAR 

WRITE (12, 7  0  0 ) Y , BETAD , RTPT , GPT 

700  FORMAT (2X,6HFREQ  = , F7 . 0 , 3X , 1 1HBETA ( DEG)  =,F7.2,5X, 
16HGAIN  =,1PE12.4, 1X,2H0R, IX, 1PE12.4 ,2HDB) 

WRITE ( 1 2 , 70 1 ) OX ,  OY 

701  FORMAT ( 2X , 6HOMEGA= , F7 . 2 , 3X , 12HATTENUATION= , F7 . 2 ) 

DBT (M) =  GPT 

OMG (M) = Y 
PHASE (M) =BETAD 
IF(NNN.GT.l)  GOTO  80 
WRITE (12,810) 

810  FORMAT ( 1H0 , 13X , 6HLENGTH , 9X , 8HDIAMETER, 7X , 8HPRESSURE , 
19X , 5HC  ADB , 10X , 3HFN  ,  7X , 7HDENSITY , 5X , 11HREYNOLDS  NO 
2) 

802  FORMAT (1H0, IX, 5HLINE  ,I2,3X,  8 ( 1PE12 . 4 , 2X)  ) 

DO  803  1=1, N 

803  WRITE (12,802) I, AD (I) , DI ( I ) , P ( I ) , CA ( I ) , FN ( I ) ,RHO(I) , 
1REY(I) 

NNN  =  2 

80  CONTINUE 

GOTO  (85,611) , IND 
85  Y=0. 

DW=  5  . 

NST  =  200 
IND  =  2 
GOTO  40 

C  THE  FOLLOWING  SECTION  IS  FOR  THE  CALCOMP  PLOTTER 
611  IF(ICAS) 502,502,610 
610  CALL  FACTOR (0.625) 

OMG (M+l) =0 . 

OMG ( M  +  2 ) =100. 

DBT (M+l) =-30  . 

DBT (M+2) =5.0 

CALL  AXIS ( 0 . , 0 . , 17HFREQUENCY  (HERTZ )  , -17 , 10 . 0 , 0  .  , 
10MG (M+l) , OMG (M+2) ) 

CALL  AXIS(0. ,0. ,1 5HGAIN  ( DECI BELS )  , 1 5 , 8  .  ,  9 0  .  , 

1DBT (M+l) , DBT (M+2) ) 

CALL  LINE (OMG , DBT , M , 1 , 0 , 4 ) 

READ (11,*) CASE (1) 

READ (11,*) NPTS , LSMB 
WRITE (12,207) 

20  7  FORMAT (1H0 , 10X , 1HN , 4X , 4HFREQ , 1  OX , 2HPS , 10X , 2HPR , 

1 10X , 4HGAIN , 10X , 5 H PHASE ) 

J  =  0 

DO  69  1=1, NPTS 
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612  READ (11,*) FREQ , PS , PR , PHT 
GPX  =  20 . * ALOGIO ( PR/ PS ) 

PHD = -PHT* FREQ* .36 

WRITE (12 ,209) I , FREQ, PS , PR,GPX , PHD 
209  FORMAT (7X,I5,5F12.5) 

DB ( I ) =GPX 

OMGX(I) =FREQ 

IF(PHD.GT.O)  GOTO  69 

J  =  J+1 

PH ( J) =PHD 

OMGP ( J) =FREQ 

JMAX= J 

69  CONTINUE 

OMGX (NPTS+1) =OMG (M+l) 

DB (NPTS+1 ) =  DBT (M+l) 

OMGX (NPTS+2) =OMG (M+2) 

DB (NPTS+2) =DBT (M+2) 

CALL  LINE ( OMGX , DB , NPTS ,1,-1, LSMB ) 

CALL  PLOT (15. 0,0, -3) 

CALL  FACTOR (0.625) 

OMG (M+l) =0  . 

OMG (M+2 ) =100. 

PHASE (M+l) =-360. 

PHASE ( M+  2 ) =45. 

CALL  AXIS ( 0 . , 0 . , 17HFREQUENCY  (HERTZ )  , -1 7 , 10  .  ,  0  .  , 
10MG ( M+ 1 ) , OMG (M+2) ) 

CALL  AXIS (0. ,0. ,21HPHASE  ANGLE  (DEGREES) , 21 , 8 .  ,  90 .  , 
1 PHASE ( M+ 1 )  , PHASE ( M+  2 )  ) 

CALL  LINE (OMG, PHASE, M, 1,0, 4) 

OMGP ( JMAX+ 1 ) =OMG (M+l) 

OMGP ( JMAX+2) =OMG (M+2) 

PH ( JMAX+1) =PHASE (M+l) 

PH (JMAX  +  2) =  PHASE (M+2) 

CALL  LINE (OMGP , PH , JMAX ,1,-1, LSMB ) 

CALL  SYMBOL (15. ,5. , .165, 5HCASE  ,0.,5) 

CALL  SYMBOL (16. ,5. , . 16  5,CASE(1)  ,0.  ,2) 

1003  CALL  PLOT ( 10 . 0 , 0 , -3 ) 

501  CALL  PLOTE 

502  STOP 
END 

Q*  **************************  * 

SUBROUTINE  RTCMP ( X , Y , A , B ) 

CALL  ANGL ( TEMP , A , B ) 

TEMP= . 5*TEMP 
Y=A*A+B*B 
X=SQRT (Y) 

X=SQRT (X) 

Y=X*SIN (TEMP) 

X=X*COS (TEMP) 

RETURN 

END 
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SUBROUTINE  CMPMP (X,Y,A1,A2,B1,B2) 

X=A1*B1-A2*B2 

Y=A1*B2+A2*B1 

RETURN 

END 

SUBROUTINE  CMPDV (Cl , C2 , A1 , A2 , B1 , B2 ) 

TEMP=B1*B1+B2*B2 

C1=A1*B1+A2*B2 

C1=C1/TEMP 

C2=B1*A2-A1*B2 

C2=C2/TEMP 

RETURN 

END 

SUBROUTINE  HSINX(ARG,X) 

A=EXP (X) 

B=EXP ( -X) 

A  =  A-B 
ARG= . 5*A 
RETURN 
END 

SUBROUTINE  HCOSX(ARG,X) 

A=EXP (X) 

B  =  EXP  ( -X ) 

A  =  A  +  B 
ARG= . 5*A 
RETURN 
END 

SUBROUTINE  ANGL (C , A , B) 

DATA  PI/3.1415926/ 

C=ABS (B/A) 

C=ATAN (C) 

IF  (A.GT.O.)  GOTO  10 
I  A=  1 
GOTO  7 
IA  =  0 

IF  (B.GT.O.)  GOTO  10 
1B  =  2 
GOTO  15 
IB  =  0 

I A= I A+ IB+ 1 

GOTO  (35,30,25,20) , IA 

C=C-PI 

GOTO  35 

C  =  -C 

GOTO  35 

C=PI-C 

RETURN 


SUBROUTINE  CALZIN (AZINI ,BZINI ,AZRNI ,BZRNI , AZIN2 , 
1BZ IN2 , AN1 , DI 1 , BTN1 ) 

TEMP=AN1*DI 1 

CALL  HCOSX (ARG1,TEMP) 

CALL  HSINX (ARG2 ,TEMP) 

TEMP=BTN1*DI 1 
CSBI 1 =COS { TEMP) 

SNBI1=SIN (TEMP) 

ZR=0  . 

CALL  CMPMP (TEM1,TEM2,AZIN2,BZIN2,ARG1,ZR) 

CALL  CMPMP (TEM3 , TEM4 , AZRN1 , BZRN1 , ARG2 , ZR) 

CALL  CMPMP ( TEM5 , TEM6 ,AZIN2,BZIN2, ARG2 ,  ZR) 

CALL  CMPMP (TEM7 ,TEM8 , AZRN1 , BZRN1 , ARG1 , ZR) 

A1=TEM1+TEM3 

B1=TEM2+TEM4 

A2=TEM5+TEM7 

B2=TEM6+TEM8 

CALL  CMPMP (TEM1 ,TEM2 ,A1 ,B1, CSBI 1 , ZR) 

CALL  CMPMP (TEM5 ,TEM6 ,A2 ,B2 , CSBI 1 , ZR) 

CALL  CMPMP (TEM7 , TEM8 , A1 , B1 , ZR,SNBI1) 

CALL  CMPMP ( TEM3 , TEM4 , A2 , B2 , ZR , SNBI 1 ) 

TEM3  =TEM3  +TEM1 
TEM4=TEM4+TEM2 
TEM7=TEM7+TEM5 
TEM8=TEM8+TEM6 

CALL  CMPDV (TEM1 , TEM2 , TEM3 , TEM4 , TEM7 , TEM8 ) 

CALL  CMPMP (AZINI ,BZ INI, TEM1 ,TEM2 ,AZRN1 ,BZRN1) 

RETURN 

END 

SUBROUTINE  ZEBRA ( C 1 , C2 , A1 , A2 , B 1 , B2 ) 

Dl  =  1 . 

D2  =  0  . 

CALL  CMPDV ( ARG1 , ARG2 ,D1,D2,A1,A2) 

CALL  CMPDV ( ARG3 , ARG4 ,D1,D2,B1,B2) 

ARG1 =ARG1 +ARG3 
ARG2  =ARG2  +ARG4 

CALL  CMPDV (Cl ,C2 ,D1 ,D2  f ARG1 , ARG2) 

RETURN 

END 


Appendix  D 


Experimental  Equipment 


Item 

Make 

Model  Serial  Numbc 

1) 

Pneumatic  Driver 
Assembly 

Bendix 

PC-DCC 

-PS 

001 

2) 

Wave  Analyzer 

Hewlett 

Packard 

302A 

018-02012 

018-01985 

3) 

Osci  1  loscope 

Tektronix 

7613 

4H4046 

4) 

Pressure  Transducer 
Dynamic 

Kr istal 

701A 

65079 

65080 

5) 

Pressure  Transducer 
Static  (0-50  psig) 

Bell  & 

Howel  1 

1000-02 

9302 

6) 

Charge  Amplifier 

Kistler 

568 

704 

1429 

7) 

Digital  Multimeter 

Hewlett 

Packard 

3466A 

3230 

8) 

DC  Power  Supply 

Hewlett 

Packard 

6205B 

4646 

9) 

Frequency  Counter 

Computer 

Measurement 

726C 

ENY0035 

10) 

Rotometer 

Brooks 

Sho-Rate 

9175-83 

11) 

Bell  Prover 

Coll ins 

Chain 

Compensated 

Gasometer 

2683 

Vita 
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